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Abstract 

The well-studied phenomenon of X-ray edge singularities is revisited 
using the sea-boson approach that has recently been placed on a rigorous 
footing. We are able to reproduce the well-known result namely, Mahan's 
power law divergences. Unlike the work of Schotte and Schotte, no lin- 
earization of the bare fermion dispersion is needed, which, by their own 
admission, is a source of some difficulty. Our approach also brings out 
some differences between the different dimensions which is not present in 
their work. Finally, our work also allows for easy generalization to poten- 
tials more realistic than the simple delta-function used commonly in the 
literature. 



1 Introduction 

The X-ray edge problem is by now extremely well understood. The work of 
Anderson ||l|,Mahan 01 and Nozieres et. al. [||, are the early pioneering works. 
Schotte and Schotte]^ were the first to apply the technique of bosonization to 
this problem. More recently, Efiler and Frahm Q have used the exact Bethe 
ansatz solution to the Hubbard model in ID, in their study of the X-ray edge 
spectra. The striking feature of this phenomenon is the power law divergence 
of the absorption close to and above a certain threshold. This is known to 
come about as a result of the filled Fermi sea being strongly perturbed by the 
emergence of the core hole, thereby a significant departure from the case when 
the hole is absent may be expected. The precise nature of this departure requires 
detailed analysis using the techniques of Many-Body theory to which we now 
turn. The sea-boson technique is so versatile that it allows for the study of this 
problem with varying degrees of sophistication and rigor. The simplest version is 
the one found most commonly discussed in the literature. It involves assuming 
that a core hole appears instantaneously when x-ray photons are absorbed, 
with a constant (time independent) interaction between the electrons in the 
conduction band and the core hole. Clearly, this is an over simplification. A 
more desirable model involves a time-dependent build up of the interaction 
between the electron and the hole that reaches a final screened equilibrium form 
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determined dynamically. This sort of build up of screening has been discussed by 
us in the past in a different context 0. In this article, (with Roman numeral I) we 
consider the simplest version as outlined in the work by Schotte and SchotteQ. 
In a future publication, we intend treating the more realistic version. 



2 The Hamiltonian 

The hamiltonian in the Fermi language may be written down as follows. We try 
to follow the notation of Schotte and SchotteQ. In what follows, Ck annihilates 
a conduction electron with momentum k and b annihilates a core electron (not 
a hole). 

^ = E^k4ck-^ 5] vik-k)cic^,bb^ + Eob^b (1) 

Being fermions, we expect the operators to obey, {ck,cjj.} = 1 and {b,b^} — 1 
and all other anticommutators involving any two of these operators are zero. 
Here Ck = k"^ /{2m) is the kinetic energy of the conduction electrons and m > 
is the band mass of the electron. We use units such that h = 1. Also 
£^0 < is the single nondegenerate core energy level, which is occupied by 
the core electron, unless excited. We assume that the intra-band interactions 
may be ignored, since the important physics is contained in the inter-band 
interaction. We assume, just as all others do, that the interaction between the 
core electron and the conduction electron is of the delta-function type which 
makes the analysis simpler. However, one of the main selling points of our 
technique is the ease with which one may study more realistic potentials in 
three dimensions as well as one and two dimensions. In the simplest case we set 
= "^o- The way Eq.( ^ is written, it is clear that vq > 0. We now proceed to 
solve for the absorption spectra using Fermi's Golden rule. 



3 Absorption Spectra from Fermi's Golden Rule 

We borrow from the work of Schotte and Schotte Q who have computed the 
absorption from Fermi's Golden rule. 




The initial state is |i„) consists of a noninteracting filled Fermi sea with n 
particles. It is the ground state of the hamiltonian below. We have set b^ = 1 
and = 0. 

H^^Y1 ^kc|cCk + Eo (3) 

k 
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The final states are the eigenstate of the hamiltonian that is shifted correspond- 
ing to the emergence of a hole in the core state and the core electron being 
promoted to the conduction band. The final state is an eigenstate of a set of 
n + 1 electrons interacting with a core hole. Here we have to set 6^6 = and 
66t 1. 

Hf=j:^Ac^-T.'^P- (4) 

k q/O 

Define, 



and, 

J^it) = {i\e' * a e-' * \i) (6) 

It can be shown quite easily that, 

/•OO 

W{lj) = 2n Ira- / e*'^* T{t) dt (7) 
I" Jo 

The problem is now well-posed and we proceed to solve it using the sea-boson 
method. In our earlier workQ, we have introduced the sea-boson method and 
made it rigorous in another recent preprint Q. The sea-displacement operators 
^k(q) are introduced and postulated to be exact bosons in these sense below. 

[Ak(q), A'(q')] = ; [Ak(q),4,(q')] = ^k.k"5q,qm;^(k-q/2)(l-n^^(k+q/2)) 

(8) 

The Fermi bilinears are related to these objects as follows( q 7^ ). 

- nj.(k) + ^Lq/2(q)^k-q/2(q) - ^I+q/2(^)^k+q/2(q) (9) 

q q 

4+q/2Ck-q/2 = A(-q) + 4(q) (10) 

In the sea-boson language, we have, 

- ^c.k(q)4(q)^k(q) +£^0 (11) 

k,q 

Hj^Y. '^k(q)4(q)^k(q) - E ^ [^k(-q) + v4[(q)] (12) 

k,q k^q^tO 

Here, u.'k(q) = k.q/m. Following Schotte and Schotte we introduce a unitary 
transformation that transforms Hi into H f 

Hf^S^HiS- e'o (13) 
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S = exp 



E 

kq 



N 



where, 



kq 



Wk(q) 



vl Ak(-q) 
iV2 c^k(q) 



(14) 



(15) 



and Ak(— q) = nj?(k — q/2)(l — np(k + q/2)). The field operator may also be 
expressed in terms of the sea-bosons as we have showed earlier p|. 



n(x)=5]e''i 



ipi'x.j = e 
1 



-i n(x) i *([pl;x) / 



(16) 



-2 i A^Co 



^t^k(q) Uk(q)+4(-q) 



(17) 
(18) 



<i>([p];x) =^e-'i- C/q ^ [Ak(-q)+4(q) 
q k 

Here ?7q is a real number whose exact form is not important. It is clear from 
the definition that a — ?/;(0)/y/po, where po — N/V is the density of conduction 
electrons. This means, 



W{u) = 2tt Im- / dt e*("+^o)t {i\B'^ {t) B{Q)\i) (19) 
Jo 



Where, 



B\t) = 

It can be shown after some algebra. 



{i\B\t) B{0)\i) = exp 



^|Tk(q)P Ak(-q) (e- --('1) * - 1 

kq 



rr , ^ ^k(q) ^q , . rr 



(20) 



(21) 



(22) 



The quantity C/q may be eliminated by making contact with free theory. Using 
the approach outlined in our earlier work|^, we multiply and divide Eq.( pT| ) 
by the noninteracting version of the correlation function, in the numerator, we 
use the expression deduced using Fermi algebra and in the denominator, use the 
form in Eq.( pi|) but with = 0- Thus we have. 



{i\B\t) B{Q)\i) = exp 



5]Kk(q)P Ak(-q) (e--'<(i)*-l 

kq 



(*|i?^(t) i?o(0)K) 
(23) 
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wk(q) Vq Y , tt2 ( '^k(q) ^ ^ 



here, 



Miraculously, the dependence on t/q cancels. Also, the free case is evaluated 
quite simply as, 

(^\Bl{t) So(O)K) = l(V(0,i)VnO,0)) = ^ - «^(k))e-^^''* (25) 

The best way to evaluate these expressions is to work with the Fourier transform 
of the exponent with the double summation. Define, 

J(w) = |tk(q)pAk(-q)<5(u; - ^) (26) 

kq 

This may be rewritten as, 

fM = ^Y. Al.(-<1)*("' - ^) (28) 

k 

It is clear at the outset that since Ak(— q) = if k.q < 0, we must have 
/(q, w) = if w < 0. Since we are interested in the t ^ oo limit, we have to 
investigate the a; — > 0+ limit of these expressions. 



/(q,a;) = ^Y.^{kF-{k' W l^-m ^)^)e{{k'' W I'^+m u:)"^ -kF)5{^-^) 

k 

- ^^''^ j dkk'^-^e{kF-{k'' + q'/A-muj)^)e{{k'' +q^l4. + muj)^ -kp) 
Jo 



Po Jo 



X / d^d 6(ij - — ) 
J m 



3.1 Computations in ID 

Here we choose = vq- This enables a comparison with the work of Schotte 
and Schotte. We may evaluate in 1D( Since fc > ), 



/ 



aild^i d[LO )= d(u> — 



m m 
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fiq,u)^0{Lu) 



2kF\q\ 



-+-^ — m uj)^) 0{{ — - — h-V+"i uj)^~kF) 
4 a 4 



Unfortunately, f{q,Lu) is not an analytic function of lu and it is not possible 
to expand in powers of to. We will have to evaluate the sums f{q,uj) and 
fiq,i-i!)/eq. For this we have to identify the regions in g € [0,cx)) for which 
the product of the theta functions is equal to unity. This is done by equating the 
arguments of the theta functions to zero and computing the roots. There are 
four roots, each pair defines an interval of integration. The intervals in question 



are q € [^(1 - ^),ff (1 + ^)] and q e [2kF - ,2fcf^ 
order for our theory to be consistently, we have to assume |g| 
ignore the second interval. 



7nuj 1 



<< kf 



Since in 
we shall 



N ^ 
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![q,uj) = uj 0{uj) 



kp 



N ^ 

1 



n( ^ 2™ 1 



Therefore, 



where. 



52 _ A- 

TT^ IT 



s_ 

n 



m vq 



(29) 
(30) 

(31) 
(32) 



Qa[(jj) = — dk diuj - —) = 
kp Jkp 2m 



A similar approach is called for when evaluating the corresponding coefficient. 

- ep) (33) 

fi.p III 

The Fourier transform may be written as, 

(34) 

Lb 

g{t) cx e^(*) g^it) (35) 
Log (\\ (36) 



1 V2to tp e""^* 
kp m it 



m 



In other words, 



foa cos 

W{uo) (X dt 

^0 



[lo + Eq — ep)t 



cx (c^ + £;;-ef^)-[25A-(V-)^l (37) 



This is in exact agreement with the result of Schotte and SchotteQ. 
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3.2 Computations in 3D 



Let us now study the screened long-range interaction. This includes implicitly 
both the genuinely long-ranged and the genuinely short-range cases as limiting 
cases. We shall soon see that the unscreened interaction leads to a divergent 
absorption for all frequencies, hence it is important to screen the long-range 
interaction. We choose the simple Thomas-Fermi screening. In 3D, we have, 

m / m Lo\ „,kp\q\ , , . 

/(q, uj) = -- e{^^ - uo)e{qma. - kl) (38) 



4vr2 po V k 

Here we have assumed tacitly, oj « ep and \q\ « kp. It is clear from the 
general definition that /(q, w) = for |g| > ikp. However, since our approach 
is vahd for small \q\ « kp, we set qmax = kp/i, following the suggestion of 
Schotte and SchotteB. Here we have, 



If we choose Vq = {^i^poe^ /f-oo){q'^ + qTF)~^i is clear that for a; — > only 
the positive part in the above sum contributes just as earlier, and we expect an 
approach to zero at the threshold rather than a divergence. 

I{u;) = ^ (40) 

where Sq = {'KaBqTF)~^ and a^^ = me^/coo- This means. 



cos 



{u + Eo- ep)t 



W{u)^ j dt '-<^{uj + E,-epYo (41) 







4 Conclusions 

We conclude by examining these results and comparing them with the results of 
Schotte and SchotteQ. In ID, we have found nearly exact agreement with their 
results. The only difference is that the exponent seems to be different somewhat. 
In our case for a fixed vq and m the exponent scales as fc^^, whereas in their case 
it seems to scale as fc^^. This is not surprising since in the Tomonaga model we 
have to choose kp oc m since we take kp^ m — > oo such that kp/m — vp < oo 
in order to make the dispersion linear. In Fig. 1 we see the various cases in 3D. 
For strong screening, we have infinite slope and for weak screening we have a 
zero slope and for optimum screening we have slope equal to unity. 
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Figure 1: Threshold X-ray Spectra 
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